Spin- Transfer Torque in Helical Spin-Density Waves 



OO 

o 
o 

(N 

Q 



O. Wessely/'^'^ B. Skubic,'^ and L. Nordstrom'^ 

^Department of Mathematics, Imperial College, London SW7 2BZ, United Kingdom 
^Department of Mathematics, City University, London ECl V OHB, United Kingdom 
^Department of Physics, 



University, Box 530, SE-75121, Uppsala, Sweden 
(Dated: December 13, 2008) 

The current driven magnetisation dynamics of a helical spin-density wave is investigated. Expres- 
sions for calculating the spin-transfer torque of real systems from first principles density functional 
theory are presented. These expressions are used for calculating the spin-transfer torque for the 
spin spirals of Er and fee Fe at two different lattice volumes. It is shown that the calculated torque 
induces a rigid rotation of the order parameter with respect to the spin spiral axis. The torque is 
found to depend on the wave vector of the spin spiral and the spin-polarisation of the Fermi surface 
states. The resulting dynamics of the spin spiral is also discussed. 
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The spin-transfer torque (STT) provides a means of ma- 
nipulating the magnetisation of a material using a cur- 
rent. The effect was proposed from theoretical consid- 
erations by Slonczewski^ and Berger^ and has attracted 
a lot of attention due to its potential use in applications 
where a magnetic state is altered by a current as opposed 
to traditional techniques involving magnetic fields. 

Most theoretical and experimental work on the STT 
concerns layered magnetic structures where the STT ap- 
pears as a consequence of angular momentum conserva- 
tion as a current traverses an interface between two re- 
gions of non-parallel magnetisation. Recently^ we have 
shown that a STT also occurs in systems with a heli- 
cal spin-density wave (SDW). This was shown by a first- 
principles calculation of the STT for a bulk rare-earth 
system (Er) in its low temperature helical spin spiral 
(SS) structure. The effect of the current induced STT 
in Er is a rigid rotation of the SS order parameter. This 
phenomenon is a bulk effect in contrast to STT in lay- 
ered systems where the STT mainly occurs close to the 
interfaces between the layers. 

In this Article we present a more general discussion 
on the current induced torque in a SS. First the effect 
is illustrated using a one-dimensional model of a SS. We 
in the subsequent section present first-principles calcula- 
tions of the STT for three real systems, Er and fee Fe 
or 7-Fe at two different lattice volumes. Previous first- 
principles calculations have shown that both Er and 7-Fe 
have helical SDWs4i^ In many other respects these two 
systems are very different. Erbium is one of several rare- 
earth elements that exhibit a helical SDW and where the 
ordering is driven by a nesting between parallel sheets of 
the Fermi surface (FS). Iron belongs to the 3d-transition 
metals where SDWs are less frequent. The fee phase of 
Fe, which is the phase where a SS order has been found 
both in experiments and theory, has only been stabilised 
at very special conditions, and it is believed that the 
helical SDW is the combined result of several parts of 
the FS. For the one-dimensional model, as well as the 
considered real systems, we find that an applied current 




FIG. 1: Local magnetisation direction of a helical spin density 
wave. The spin density wave is characterised by the SS wave 
vector q, the cone angle 9, and the rotation angle = q ■ r. 



induces a STT generating a rigid rotation of the spiral 
order. With this theoretical prediction, new types of po- 
tential applications of the STT can be imagined, such as 
current driven oscillators with tunable frequencies. The 
current driven spin dynamics for such a device is investi- 
gated in section (|VIip 



II. THEORY 

For a system with a helical SDW, the direction of the 
local magnetisation rotates around the SS quantisation 
axis as one moves in the direction of the SS wave vector. 
Without loss of generality, we here consider the rotation 
axis to be parallel to the SS wave vector, except when 
explicitly stated, and will refer to it as the SS axis. Al- 
though we here consider helical SDWs, the results are 
valid for cycloidal SDWs as well. A SS is characterised 
by its SS wave vector q and its cone angle 9, which is 
the angle between the SS axis and the local magnetisa- 
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tion, see Fig. [TJ In this section we first consider a one- 
dimensional model of a SS, where the spin currents and 
STT are calculated analytically. We then discuss how 
accurate first-principles calculations of the STT can be 
performed. 



A. The STT in a ID model of a spin spiral 

As a simple one-dimensional single band model of a 
SS system, we consider independent particles in a spin- 
dependent potential. The spin-dependent potential is 
chosen such that the ground state of the system has a 
helical magnetisation M = Mm where 



m — (sin 9 cos qz, sin 9 sin qz, cos 



(1) 



Here, the SS wave vector q — qz and 9 is the cone angle. 
The single particle Schrodinger or Kohn Sham equation 
for such a model system takes the form 



between zq and zq + dz is given by the spin flux into the 
region 



dJ dQ 



dz. 



(6) 



Both the spin- flux 83/ dt and the current j depend on 
the change in occupation of the states at the Fermi level 
caused by the applied electric field. These two quantities 
are linked by the expression. 



dJ 



(7) 



which defines the torque current tensor C. The torque 
current tensor is in three dimensions a matrix, but be- 
comes a vector in a one-dimensional system, where also 
the current vector is reduced to a scalar. For our model 
system, the torque current tensor can be calculated by 
inserting the Fermi Bloch states, ipkp and ip-kp from 
Eq. ^ into Eq. (gl)-®, resuhing in, 



2me dz^ 



■ [sin 6 (cos qz ax + sin qz ay) 



-cos 9 



4^ 



= 0, 



where U is the exchange splitting, rrie is the electron 
mass, fi is the Planck constant ax^ ay, a^ are the Pauli 
matrices and e is the one particle energy. This equation 
can be solved analytically^ and the eigenfunctions are the 
so called generalised Bloch states. 



■0^ 



i'nk{zA) 
i'nk{zA) 



cos (0fe/2) e-^9-/2 , 

lqz/2 ] ■ l^J 



sin {9k 1 2) e 



For the above wave function the direction of the spin s 
rotates around the SS axis, with rotation angle (j) ~ qz, 
at an angle 9k with respect to the SS axis, 

s(s) = (sin 9k cos qz, sin 9k sin qz, cos 9k). (3) 

Note that the single electron cone angle 9k and the cone 
angle 9 of the SS, see Fig. [Tl in general are different. In 
one dimension, the charge current j, induced by applying 
an electric field, is a scalar given by 



-Sf—ReUl 
rup V 



■ ^ / 
oz 



• ^ / 



(4) 



where 5f is the change in occupation at the Fermi level 
due to the electric field, fej? is the Fermi wave vector and 
e is the elementary charge. In a similar way, the spin 
current Q due to a weak electric field is a vector given 
by 



Qi{z) = -Sf -^ReUl^i-^aipki. 



(5) 

From the spin current can the rate of change of angular 
momentum within an infinitesimal region of the system 
be calculated. The rate of change or torque on the region 



9Q 

dz 



dz/j{za) = a{9,kF) 



-q sin (gzo) 
q cos (gzo) 




dz. 



(8) 

where the first factor, a{9, kp), depends on the spiral cone 
angle 9 and the Fermi wave vector kp. For planar spin 
spirals, where 9 — 7r/2, the factor reduces to a simple 
function of the single electron cone angle 9k of the states 
at the Fermi level given by 



a(n/2,kp) — — — ; — . 

^ ' ' 2e 1 - (g/2fc)cos 9kp 



(9) 



A comparison between Eq. ^ and Eq. ^ shows that the 
current induced torque, in our one-dimensional model, is 
perpendicular to the magnetisation direction and lies in 
the plane of rotation of the SS, thus causing the SS to 
rotate with respect to the SS axis. It can also be seen 
from Eq. ([5]) that the magnitude of the STT is governed 
by the length of the SS wave vector q and the factor 
a{9, kp). 



B. The STT for a real system 

The STT for a real system can be calculated by general- 
ising the expressions for the charge and spin currents in 
one dimension to three dimensions. In three dimensions 
can the currents induced by an external electric field, to 
linear order in the field, be calculated from a FS integral. 
For a multiband system must the contribution from all 
bands crossing the Fermi level be taken into account. In 
the previous section was a model system and the STT on 
an infinitesimal region at an arbitrary point of the sys- 
tem considered. For a real material, where most of the 
angular momentum is localised on the atoms is the time 
evolution of the atomic moment governed by the STT on 
the corresponding atom. The torque on an atom is given 
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by the spin-flux into a sphere surrounding the atom. For 
a single electron state with band index n and wave vector 
k, the flux into a sphere of radius R is given by, 



Q„k ■ rR sin 6 



dJr 



dt 



(10) 



In three dimensions is the one electron spin current tensor 



Q„k(r) = Re s ® v0„k(r)| , 



(11) 



where s — {h/2)a^ and v — {—ih/rrie)'^ ■ The total spin 
current is the sum of the spin currents from all occupied 
states. A external electric field E induces a non equi- 
librium spin current by changing the occupation of the 
states Sf at the FS. The change in occupation number 
df can be calculated using the relaxation time approxi- 
mation and semiclassical Boltzmann theory^ 



<5/(k) - -<5(e„k 



, re 

fF)^ VkEnk • E, 



(12) 



where t is the electron relaxation time and e„k is the 
band energy. The total non equilibrium torque on an 
atom can be calculated from a FS integral, 



d3 



-VcTe 



dt (27r)3?i 



E 



^(Vke„k-E)^:^, (13) 
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where Vc is the volume of the system and a summa- 
tion is performed over all bands crossing the Fermi level. 
The above equation defines a linear relation between the 
torque and the external field, 



93 V- . ^ 



(14) 



A similar relation can be obtained for the charge current 
density and the electric field. 



(15) 
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/ VkEnk (Vkfnk • E) 
JFS 



nk 
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where B„ is the conductivity tensor for band n. By 
combining Eqs. (fT4|) -([T5 |) is a linear relation between the 
torque and the current density obtained, where the un- 
known electron relaxation time r has been cancelled, 



d3_ 

dt 



(EAn)(EM^'j = Cj. 



(16) 



The above equation defines the spin current tensor C in 
3 dimensions. 



C. The STT from the augmented plane wave 
(APW) method 

The wave functions used in the previous section to cal- 
culate the torque current tensor C can be obtained from 
first principles electronic structure calculations using spin 
density functional theory. Several electronic structure 
methods use basis sets where space is divided into muffin- 
tin spheres surrounding the atoms and an interstitial re- 
gion. This division is introduced in order to simplify 
the calculation of the intra and interatomic behaviour of 
the wave functions. When calculating the STT, a nat- 
ural choice is to use the atomic augmentation (muffin- 
tin) sphere as the surface for evaluating the STT on the 
atoms. We now illustrate more in detail how the calcula- 
tion may be carried out within the augmented plane wave 
method (APW). The APW expansion can be written as 
a sum of plane waves 



i(G+k-q/2)r 



Mr) - 2^ \ ^^^^^g.(G+k+q/2)r 



(17) 



where G are the reciprocal lattice vectors. The plane 
wave coefficients a„k and 6„k are obtained from the first 
principles calculation. The spin-flux into a sphere with 
radius R centred at an atom at site r„j is for plane waves 
given by the expression 



^0 



Qnk • f i? sin 9 d(j)d9 



hR^ 



Re ^ -Utt 



G,G' 



a^sa a,*k,Gank.G'e-*(°-°')"--ji(|G - G'\R) (G' + k - q/2).(G^G') 

ats/3 a:k,G^.k.G'e-'(°-°'-'»'"-'"Ji(|G - G' - q|i?) (G' + k + q/2).(G -'G^ - q) 

/Jtsa 6:k,Ga«k.G'e-''(°-°'+'i)--'"ji(|G - G' + q|i?) (G' + k - q/2).(G -'G^ + q) 



/3ts/3 6:k,Gfenk,G'e-^(°-°'>"ji(|G - G'\R) (G' + k + q/2).(G^G' 

I 



(18) 



where ji is the first spherical Bessel function and a, (3 pression evaluated it is straight forward to calculate the 
are the up and down spinors respectively. With this ex- 
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STT using Eq. ^ and Eq. 



III. MATERIALS WITH SPIN SPIRAL 
MAGNETIC ORDER 

Long ranged magnetic order in form of helical SDWs ex- 
ists in several types of materials. Maybe the most fa- 
mous is the heavy rare earth elements. They have similar 
valence electron structure, resulting in similar chemical 
structure, and all the tri-valent elements order in the hep 
lattice structure with similar lattice volumes. The main 
difference along the series lies in the filling and magnetic 
moment of the localised 4/-electron shell. The FS of Er 
which is typical for the series has a strong nesting feature, 
i.e. two large parallel sheets of the FS. The formation of 
a SDW with a wave vector equal to the nesting vector 
allows for hybridisation removing these parts of the FS 
opening up a gap at the Fermi level lowering the total 
energy of the system. Anisotropics in the system deter- 
mine whether the type of SDW preferred in the system, 
is a helical SDW, conical SDW or a longitudinal SDW. 

Helical SDWs also exist for various 3d-transition com- 
pounds, however among the elemental 3d-transition met- 
als the magnetic structure are less exotic at ambient con- 
ditions. It was previously shown that if certain general 
criteria are met, helical SDWs are formed in the transi- 
tion metal series. This means that at conditions slightly 
different from ambient, spin spirals occur naturally in- 
dicating that this type of magnetic order is less exotic 
than one could believe. An example of such a system 
is 7-Fe the fee phase of iron. Experimentally 7-Fe has 
been stabilised as precipitates in a Cu-matrix. The ex- 
perimental magnetic ground state structure was found to 
be a SS with q = 27r/a(0.1, 0, 1) where a =6.76 a.u. Sev- 
eral theoretical studies^'^^.'-^ have mapped out the ground 
state magnetic structure of 7-Fe which was found to be 
very volume sensitive. In contrast to the rare-earth sys- 
tems where the SS state is driven by nesting between two 
parallel sheets of the FS, the SS state of 7-Fe is not pro- 
moted by a single FS nesting vector. Instead there is a 
net energy gain from many parts of the FS given by the 
hybridisation by a SS vector. 



IV. ELECTRONIC STRUCTURE 
CALCULATION 

All the materials specific parameters for Er and 7-Fe used 
for the calculation of the torque current tensor C were 
obtained from first principles density functional theory. 
The calculations were made using the full-potential aug- 
mented plane wave plus local orbitals (APW+lo) method 
as described in Ref. For Er, the calculations were per- 
formed in the same manner as in Ref. M and for 7-Fe 
the calculations were performed in the same way as in 
Ref. ii. The local spin density approximation (LSDA) 
as parametrised by von Barth and Hedin was used with- 
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FIG. 2: The top panel shows the total energy of Er per hep 
unit cell for SS wave vectors from F to ^. The middle panel 
shows the total energy per atom of 7-Fe for SS vectors from 
F to X in the larger lattice volume. The bottom panel shows 
the total energy per atom of 7-Fe for SS vectors from X to 
W in the smaller lattice volume. 



out any shape-approximation to the non-coUinear mag- 
netisation, i.e. charge and magnetisation densities as 
well as their conjugates potentials are allowed to vary 
freely in space both regarding magnitude and direction. 
For Er, a set of 12x12x8 k-points was used for con- 
verging the electron density and for Fe we used a set of 
20x20x20 k-points. The SS was treated using the gener- 
alised Bloch theorem. The Er calculation was performed 
with the 4/-electrons treated as spin polarised core elec- 
trons and we used the experimental Er hep lattice param- 
eters (a = 6.73 a.u. and c = 10.56 a.u.). For Er, we cal- 
culate a large number of SS wave vectors, q = 7r/c(0, 0, q) 
in the first Brillouin zone of the hep lattice, along the out 
of plane direction between F and A. The total energy for 
these calculations are presented in Fig. [21 There is an 
energy minimum for q = 0.40 which corresponds to the 
ground state SS structure of Er. For Fe we perform cal- 
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culations at two lattice volumes. At a lattice constant 
of a = 6.82 a.u., corresponding to the lattice constant of 
Cu, we do a series of calculations for qi = 27r/a(0, 0, g), 
between F and W in the first Brillouin zone, and an en- 
ergy minimum is found for q = 0.59 (see middle panel in 
Fig. [2]). At a reduced lattice constant of a = 6.66 a.u. we 
do a series of calculations for q2 = 7T/a{w, 0, 2), between 
W and X, where there is a weak local energy minimum 
'dt w = 0.48. This latter minimum corresponds to a spi- 
ral wave vector which is close to the experimental SS of 
7-Fe. Note that q2 is non-parallel to the SS rotation axis 
which is in the z-dircction for all three calculations. 



V. RESULTS 

In this section we present the results of the STT cal- 
culations for the three SS systems treated in the previ- 
ous section. The one electron spin current tensor Q„k 
was calculated using Eq. (fT5)) and the Kohn-Sham eigen- 
functions obtained from the first principles calculation. 
The first Brillouin zone was covered by a 41 x 41 x 41 
k-point mesh, in order to accurately evaluate the FS inte- 
grals in Eq. and Eq. ([TS]). The conductance tensor, 
B = J2m. ^"i, was as expected, found to be diagonal for 
all three considered systems. 

The torque current tensor C was calculated for Er hav- 
ing a SS with wave vector q — n/c{0 0.4). We found 
that for an Er atom situated at a site with magnetisation 
direction (10 0), the spin current tensor 



0.0 0.0 
0.0 0.0 
0.0 0.0 



From the structure of the C tensor, it can be seen that for 
a current in the (001) direction, a torque is induced in 
the (010) direction. Since the torque is perpendicular to 
the local magnetisation and lies in the SS rotation plane, 
we conclude that a current along the SS axis causes a 
rigid rotation of the SS around the SS axis. 

For 7-Fe with the larger volume in the qi = 
27r/a(0, 0, 0.59) spiral state, the torque current ten- 
sor, evaluated for an atom with magnetisation direction 
(100), 




0.0 0.0 0.0 
C = - I 0.0 0.0 -1.4 
^ 0.0 0.0 0.0 



Also in this case we conclude that a current along the 
SS axis causes a rigid rotation of the SS. In this system, 
the torque, for the same current density, is found to be 
nearly three times larger than for Er. We will return to 
a discussion on the magnitude of the torque later in this 
section. 

Finally we calculated the STT in 7-Fe with reduced 
volume, where the SS wave vector q2 = 7r/a(0.48, 0, 2) is 



a: Er 










n 


(?i/e)(A„)23 U 


leV] 


(RVe')(B„)33 


[eV/A] 


1 


0.01 




0.006 




2 


0.01 




0.01 




3 


-0.04 




0.03 




4 


-0.008 




0.007 




b: Fe qi 










n 


(?i/e)(A„)23 


[AeV] 


(aVe2)(B„)33 


[eV/A] 


1 


-0.002 




0.002 




2 


0.01 




0.02 




3 


-0.1 




0.07 




4 


-0.005 




0.002 




c: Fe q2 










n 


(7i/e)(A,02i 


[AeV] 




[eV/A] 


1 


0.005 




0.006 




2 


0.02 




0.01 




3 


0.02 




0.02 




4 


-0.05 




0.05 




5 


0.09 




0.08 





TABLE I: The elements of the A„ and B„ tensors that mainly 
contribute to the C tensor, see Eq. (I16|l 



non parallel to the spin rotation axis z. We found that 
the torque current tensor for an atom with magnetisation 
in the (10 0) direction is given by. 



C = ^ 



^ 0.0 0.0 0.0 

0.6 0.0 0.0 I [A^' 
^ 0.0 0.0 0.0 



The structure of the torque current tensor for this system 
differs from the previous two systems. Here, a current in 
the (10 0) direction is required to produce a rotation of 
the SS. The STT from a current in the (001) direction 
vanish since 7-Fe with q2 = 7r/a(0.48, 0, 2) is antiferro- 
magnetic in the (001) direction. 



VI. DISCUSSION 

In general, the STT depends on a balance between the 
size of the SS wave vector and the spin polarisation of 
the conduction electrons at the FS. 

Some insight into the origin of the STT on the single 
electron level can be obtained by considering the semi- 
classical expression for the spin current^i 



Q(z) = h{z)l. 
2 e 



(19) 



If the semi-classcal spin current is combined with Eq. ([S]) 
one obtains the following expression for the torque cur- 
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rent tensor. 



VII. DYNAMICS 



C = 



h ds{z) 



2e dz 



dz — — sin 

zo 2e 



-q sin (gzo) 
q cos (qzo) 

V 







(20) 

The above equation gives the same result as the quan- 
tum mechanical calculation of C in section pi A[) , if the 
denominator of the factor a in Eq. ^ is equal to one, 
which is approximately true if (g/2fc)cos 9kp is small. 
Semi-classically, the maximum spin transferred per elec- 
tron, for a given value of q, to an atomic layer is given 
by (h/2)ql where I is the inter layer distance. This maxi- 
mum can be achieved if the single electron cone angle 9kp 
is equal to tt/2. For Er and 7-Fe with the larger volume 
is the semi-classical maximum of the spin-transfer, with 
the given value of q, equal to 0.3?i and 0.9?i, respectively. 
A more accurate value for the average spin transferred 
per electron to an atom is obtained by multiplying the 
torque current tensor with e/A where A is the area per 
atom in the direction of the current. For Er, where the 
area per atom in the (001) direction A = a^\/3/2, the 
spin-transferred per electron is 0.057i. For 7-Fe, the area 
per atom in the (001) direction is and the spin- 

transferred per electron is 0.2h, for the larger volume. 
Hence, the spin-transfer torques for Er and 7-Fe (larger 
volume) obtained from the first principles calculations 
are only 16% and 22%, respectively, of the their semi- 
classical maximum. This is partially due to the fact that 
the cone angle 6kp of the conduction electrons is less than 
tt/2. 

Another reason for the reduced STT is that contribu- 
tions from different bands partly cancel each other. The 
contributions from the individual bands to the STT are 
given by the A„ and B„ tensors shown in Table For 
Er there are four bands that cross the FS and the domi- 
nant element in the A„ tensors is their (A„)23 elements. 
As shown in Table (|TK) , the STT from band three is four 
times larger and in the opposite direction to the STT 
from band one and two. For 7-Fe (larger volume) there 
are four bands crossing the Fermi level, where the dom- 
inating contribution to the STT comes from the third 
band, see Table ([I)d). For 7-Fe with the smaller volume 
there are five bands crossing the Fermi level and the dom- 
inant element in the A„ tensors is their (A„)2i elements. 
As shown in table (|TJ:), the largest contribution to the 
STT is here coming from the fifth FS. Finally we would 
like to note that semi-classically the maximal spintrans- 
fer per conduction electron is h/2, which only can be 
obtained if both ql = 7r/2 and 9kp = Tr/2. An other 
route to high spin transfer is to try to reduce the denom- 
inator in the expression for the a factor in Eq. ^ . The a 
factor becomes large for small Fermi surfaces where q/k 
is large. 



We will for our SS systems represent the uniform mag- 
netisation of an atomic plane perpendicular to the SS 
axis by a unit vector m in the direction of the magneti- 
sation. The time evolution of m is phenomenologically 
described by the LLG equationi^ii^. 



dm 
'dT 



dm 
~dt' 



(21) 



The first term T is the total torque exerted on the layer 
and the second term is the Gilbert damping term with 
damping constant a. The systems under consideration 
are all having a helical SDW with an easy-plane or hard- 
axis anisotropy. The hard-axis is directed along the SS 
axis and helps maintain the planar structure of the SS. 
For such systems is the total torque T a sum of a current 
induced STT and a torque resulting from the hard-axis 
anisotropy. Both of these torques are perpendicular to 
the magnetisation ih and the SS axis which here is con- 
sidered to be in the z direction. Eq. (PT|) can therefore 
for this system be written as 



where 



dm . ^ ^ ^ dm 

= Am X z — am x , 

dt dt 



Mj, = ^-5 + 6 m • z. 

sm a 



(22) 



(23) 



The first term in the above equation is due to the STT 
and the factor a^Ojkp) was introduced in Eq. The 
b in the second term is the strength of the hard-axis 
anisotropy. Eq. (|22p is in spherical coordinates, see 
Fig. [1] given by. 



dt 



{l + a^)^~Aij,( 



and 



d6 

_(l + a2) =-aA{j,e) sin 
dt 



(24) 
(25) 



Eq. (|24p describes a rigid rotation of the SS due to the 
current induced torque where the angular velocity, as a 
function of the current j, is given by —A{j,9)/{1 + a^). 
It can be seen from Eq. that the effect of the Gilbert 
damping is to change the cone angle of the SS. Eventu- 
ally, the spiral will reach a steady state with a cone angle 
6*0, given by 



q j 0(6*0, kp) 



-b cos 61 



sm 







(26) 



when the current induced torque and the anisotropy 
torque balance each other. At this point in time, also 
the rotational motion of the magnetic moments with re- 
spect to the SS axis stops. The polar angle (p will thus, 
when a current j is passed along the SS axis, change by 



A0 = - 



1 



(27) 
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VIII. SUMMARY AND CONCLUSIONS 



In ( tan — ^ — In ( tan — 

4 y \ 2 



(28) 



before the spiral reach the steady state with cone angle 
6q. The SS will returns to a planar state if the current 
is switched off after the steady state is reached. During 
the reversal, the spiral rotates in the opposite direction 
performing the same number of rotations as required to 
reach the steady state. 

The current induced torque has, for slowly vary- 
ing magnetisation, been described by introducing an 
adiabatic and a non-adiabatic term in the LLG 
equationi^ii^ii^iii. The adiabatic term, which for a 
SS whose axis is along z takes the form drii/dz, has 
in Eq. ((22|) been replased by the STT component of 
A(j,9). A non-adiabatic torque term proportional to 
rh X {dm/dz) could also be introduced in Eq. ([22|) . The 
main effect of such a torque would be a modification of 
Eq. (|25p , which determines the dynamics of the cone an- 
gle 9. A non-adiabatic torque would add a term to the 
right hand side of Eq. ([^5]) which could decrease or even 
change the sign of dO/dt. 

The effect of in plane anisotropy is discussed in Ref. 3, 
where it is suggested that the current needs to overcome 
a certain critical current before the spiral start to rotate. 



We have shown in general, that the a current through 
a SS induces a rotation of the spins, where the angular 
velocity depends on the magnitude of the current. The 
size of the STT has been analysed in terms of the spin- 
transfer per conduction electron and we conclude that 
the total STT depends on, the spin-polarisation of the 
electrons, the SS wave vector and how contributions from 
different bands coincide. The dynamics of the SS has 
been calculated including the effects of anisotropies and 
damping. For a system with a hard-axis anisotropy along 
the SS axis, the damping leads to a steady state where 
the rotation eventually stops. 
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